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ABSTRACT 


Three approximate methods of second-order analysis are 
compared to an "exact'' second-order analysis. They are the Fey method, 
the moment magnifier method and the amplified lateral load method. 

In general, the Fey method and the moment magnifier method 
are shown to give good approximations with a slight degree of 
conservatism. The amplified lateral load method is shown to be very 
conservative in the more flexible frames studied. 

Column design following a second-order analysis is also 
presented. It is shown that estimating the maximum design moment 
within a column is underestimated when the calculations are based on 
the braced frame effective length and that an effective length equal 


to the actual length should be used. 
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NOTATION 


cross-sectional area of bracing member 


‘modulus of elasticity 


storey or building height 

shear in a given storey 

service storey shear 

ultimate storey shear 

moment of inertia 

H/A, = lateral stiffness 

effective length factor 

length of column, center to center of joints 

length of brace 

first-order moment 

total second-order moment = first-order moment plus P-A 
effect 

larger end moment in a column, always positive 
smaller end moment in a column, positive if column is 
bent in single curvature 

axial load in a column 

critical load of a column in a storey or frame in lateral 
buckling 

T°EI/L? = Euler buckling load 

ultimate load in a column 

EPA, /Hh = stability index 

hP/EI, «7 = Tey E.) 


slope of brace in Figure 2.4(b) 
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first-order deflection 
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additional deflection due to P-A effect 
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CHAPTER I 


INTRODUCTION 


1.1 What Problem Is 

Design of a reinforced concrete frame is preceeded by a 
structural analysis to determine member end forces and moments. The 
analysis is usually effected by a computer program based on the direct 
stiffness method. The majority of such programs assume linear elastic 
behavior of all members and are first-order programs, that is, equilibrium 
is formulated on the undeformed structure. While the assumption of 
elastic behavior is reasonable for a structure under working loads, 
the neglect of the effects of deformations represents a greater 
restriction, especially in flexible structures subjected to lateral 
and vertical loads. 

An analysis which formulates equilibrium on the deformed 
structure is called a second-order analysis. If the lateral deflections 
of the isolated storey shown in Figure 1.1 are ignored equilibrium 


could be written as: 
2M = M, + M, +M +M = Hh a Ae 


Had equilibrium been formulated on the deformed storey (second-order 


solution), Equation 1.1 would have appeared as: 


sM = M, = M, so M, + M, = Hh + PA. Nia: 


The first-order solution ignores the term PA. It should also be 
noted that the deflections computed in a first-order analysis, Al» are 


less than those of a second-order analysis, A,. 
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FIGURE 1.1 DEFORMED STOREY UNDER COMBINED LOADING 
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1.2 ‘scope of Thesis 


Due to the cause and effect nature of the P-A phenomenon, 
the determination of the second-order deflections is frequently done 
in an iterative computation which can be tedious. 

The intent of this thesis is to review and compare several 
approximate procedures which obviate the need to perform successive 
iterations. Design rules based on these analyses are also verified. 

Chapter two presents a review of frame stability and second- 
order effects in frames. A series of ree frames used to compare 
the analyses are described in Chapter three. Chapter four discusses 
the approximate methods of analysis and compares the results of 
these analyses to those from an "exact" analysis. Chapter five considers 
the column design procedures which are necessary once a second-order 


analysis has been carried out. 
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CHAPTER 2 


BACKGROUND INFORMATION 


2-1 Introduction 
This chapter presents a brief review of frame stability by 
considering slenderness effects in buildings. A number of methods of 


second-order analysis are PCE ed a eos oe 


2.2 Euane stability 


Ziz.i1 Frame Response to Horizontal and Vertical Loads 


The material in this section is based on Reference 10. 


Similar analyses have been presented by eeeens Ge rey 


(18) and pai 


Goldberg 
Figure 2.1(a) represents a storey in a frame subjected to 


combined loading. Let A denote the relative deflection between the 


top and bottom floors of the given storey as a result of lateral forces 


only. The corresponding storey stiffness may then be defined as K, = H/A 


x 


in which H is the storey shear. 
Application of the vertical loads will tend to increase the 
relative lateral deflection by an amount Aa The total relative 


deflection, A is then the sum of A and A. The total storey moment 


Oe 


(sum of top and bottom end moments in all the columns) will be 


M, = Hh + PA, Para k 


in which h is the storey height and 2P is the sum of applied vertical 
loads down to the storey considered. If it is assumed that the moment 


diagrams in the columns resulting from the action of the vertical loads 
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FIGURE 2.1(a) STOREY UNDER COMBINED LOAD 
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FIGURE 2.1(b) DECOMPOSITION OF VERTICAL LOAD 
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have the same shape as those produced by the lateral forces, which 
will be essentially true if the columns are stiffer than the beams. 


The storey moments can be taken equal to hK Ay. Equation 2.1 then 


Q 
gives 
H 
A, = ‘att =P Zee 
he h 
and, 
ye 
me Sone qaqht see 2.3 
2. K. - oP 
Q h 


Hence, slenderness effects are taken into account by introducing 
a magnification factor given by Equation 2.3 that multiplies the nominal 
storey moment Hh. Reference 2 proposes that all stresses and deformations 
resulting from the sole action of the lateral forces be multiplied by 
the same factor. 

Approximate proportionality of the bending moment diagrams 
in columns, as a result of lateral and vertical forces, requires that 
the columns in every storey deflect as nearly straight lines and that 
A, be nearly proportional to A in all the storeys, so that the storey 
shear stiffness remains practically constant. 

Considering again the deformed storey of Fig. 2.1(a), if 
axial deformation and higher order terms are neglected, all joints 
in a floor will remain in a horizontal plane as the frame deflects. 
The storey drift is defined as p = A/h and is the same for all the 
columns in a given storey. The load P is now replaced by its horizontal 
and inclined components as in Fig. 2.1(b). 

The horizontal component is equal to Pw while the second is 


essentially equal to P for smail drifts and columns deflecting as straight 
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lines. Consequently, the combined action of vertical loads and an 
applied storey shear, H, on the deflected frame is equivalent to the 
action of the original axial loads on the columns, combined with an 
increased storey shear, H + Py, as salen by Equation 2.1. 

If the beams are rigid in comparison with the columns, the 
columns will not deflect as straight lines. Consider the case of 
infinitely rigid beams and the coordinate system of Figure 2.2. 

Let A(y) denote the deflections of a column, as a result of the storey 
shear alone. For a small axial load P, the bending moment increases 
by PA. The deflection of the top relative to the bottom, due to P, 

is found to be PAh*/9.87 EL assuming the column deflects as a sine 
wave. Therefore, the increased lateral deflection due to the vertical 


toad is: 


a Ph* 
A, a la, u A, OF 87 EL 2.4 


But the lateral stiffness Ko due to lateral loads only is Ko =) POU Mn 


where the sum extends to all the columns in a storey considered. From 


Equation 2.4 then, 


Q Q 
H 
A, ij a 1299" 2 2.6 
Q h 
and, 
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The term 1.22 will vary from 1.0 to 1.22 as the P-A moments 


vary from a straight line distribution to a sinusoidal distribution. 
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FIGURE 2.2 ISOLATED COLUMN WITH RIGID BEAMS 
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EDLs ss a function of i and of the relative stittnesses of the 


columns and beams. 
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The tilting effect of foundation rotation and of column 
shortening and lengthening in an unsymmetrical building may be taken 
into account by including these terms in A,. 

Zeoee rane,buckling 

The sway buckling load of a storey can be obtained by solving 

Equations 2.8 and 2.9 with H set equal to zero. 


For a non-trivial solution aLP/h must equal K, in Equation 2.8. 


Me 


Therefore, the critical load for the storey under consideration is 


xP ot Meare oes ZeL0 


since Ky = H/A, this may be rewritten as: 
ae = oe 2.11 


where @ has been taken equal to 1.0. 


2.3 Review of Second-Order Analyses 


Z,30b, Pot Treration 
In section 2.2.1 it was shown that the total storey moment 
for the frame shown in Fig. 2.1(a) was given by Equation 2.1. The term 
"P-A effect" refers to the extra bending moments which are developed 
throughout the structure when the vertical forces, XP, are displaced 
laterally through a deflection A. These extra moments are not 


determined by a conventional first-order structural analysis, since 
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such analyses formulate equilibrium on the undeformed structure. If, 
however, the shear in the storey was artificially increased by an 
amount EPA, /h, the sum of the column end moments in the storey as computed 
by a first-order analysis would be ours 

These principles may be readily extended to a more complex 
structure. Initially the structure is analyzed under the given loading 


to obtain the first-order lateral deflections, denoted as A.» A etc. 


ie 
in Figure 2.3. The artificial storey shears, which would produce 
column end moments equivalent to those caused by the vertical loads, 
are then computed: 
XP. 
' = _ = 
Me De one ee eee 
e: 
The storey shears computed by Equation 2.12 may be represented by 
additional sway forces H' which will produce moments and forces throughout 
the structure which will simulate the P-A effect. The artificial sway 


forces due to the vertical loads, Has are then computed as the 


difference between the additional storey shears at each level: 


Hei= V3 - V! peaks: 
3 i=l ab 


The sway forces, His are added to the applied lateral loads, and the 
structure re-analyzed. When the A. values at the end of a cycle are 
nearly equal to those of the previous cycle, the method has converged, 
and the resulting forces and moments now include the P-A effect. It 


should be noted that the sway forces can be either positive or negative. 


(1) 


The application of P-A analyses is discussed by Adams ; 


Co) CED) 


Wood, Beaulieu and Adams , MacGregor and many others. An 


example of the application of such an analysis to a tall steel building 


is given by Springfield and ene 
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BICURE 2.3 SWAY FORCES DUE TO VERTICAL LOADS 
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2.3.2 Simplified Second-Order Frame Analysis 


Reference 3 shows that a direct solution of the second- 
order deflections and moments can be obtained using a standard first- 
order program by inserting a fictitious diagonal brace of negative 
area in each storey as shown by the dashed lines in Fig. 2.4(a). 
The brace area may be obtained by developing the second-order stiffness 
matrix for the deformed column of Fig. 2.5. Static equilibrium of the 


column gives: 


ae M. + M : P(A, - A.) me Ae 
£ i ie b ; 


Substituting the slope deflection equations for M. and M, into 
Equation 2.14 yields the matrix of stiffness influence coefficients 


given in Equatton®* 2715. 


M, 4EL/L 2EL/L =5E T/L 6EL/L- 6. 

SLAG 4EL/L GEL le 6EI/L? 0. 

oe K—GEU/Uae orl Gs Slebl ile ee pew olny i PL A. 

Ee GEl/D we ORL lo tH 12ERL/ Leb) ee Ey = PL ae 
2.15 


A first-order analysis omits the terms involving P/L. However, 
if the structure contained bracing members as shown in Fig. 2.4(a), 
the program would generate a stiffness matrix corresponding to the 


degrees of freedom shown in Fig. 2.4(b), that is: 


F cos*R -~cos* A 
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FIGURE 2.4(a) 
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FRAME WITH NEGATIVE BRACING 
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FIGURE 2.4(b) NEGATIVE BRACING MEMBER 
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The two by two matrix would be superimposed in the same position as the 
P/L terms occupy in the second-order matrix. Considering the brace 
of the Figure 2.4(b) the required brace area is obtained by equating 


(AE/L, cos*fi) and (-P/L): 


-P Os 


L Ecos“ 


Zed 


The value of P in Equation 2.17 represents the sum of the axial loads 
in the columns of a storey. Because the brace has a negative area, it 
increases the flexibility of the frame, to account for the P-A effect. 
A first-order analysis of a structure containing these 
negative bracing members will compute forces and moments within the 
structure which correspond to the second-order A The axial force 
in the brace causes slight errors in column shears and axial loads 
which can be corrected by statics. These errors can be kept to a 


minimum by using long flat bracing members. 


2.3.5 -jsecond-Order! Kinite: Element Analysis 


K3 Ne Ses ie 


has described one type of solution of the 
second-order moments and forces by means of the finite element method. 

If a frame is visualized as an assemblage of elements inter- 
connected at their ends, called nodes, the equilibrium configuration 
of the complete structure can be expressed in terms of the nodal 


displacements. The force-displacement relationship for the element of 


Figure 2.6 can be written as: 


[k] {w} = {P} 2.18 
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FIGURE 2.6 


ELEMENT IN LOCAL COORDINATE SYSTEM 
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The second-order element stiffness matrix is comprised of two parts, 
a first-order stiffness matrix, [kK], and a non-linear geometrical 
stiffness matrix, EAR The matrix [kK] is unique to the second-order 
solution and represents the action of axial loads on a deformed 
element. Both matrices are shown in Figure 2.7. 

For an element in a global system, at an angle 9 to the local 


system, the global element stiffness matrix would be given as: 


me] = (RIT [Kd] ie] Pe) ero 


The element end forces and displacements are transformed as: 


{Pp} 


[R] {P} 


Il 


{w} [R] {w} 2.20 


The matrix [R] is shown in Figure 2.8. 

The global element stiffness matrix, transformed by Equation 
Poets Snowy in ficure..2.9. 

Similar to the force-displacement relationship for the element 
the force-displacement relationship for the complete structure, or the 


complete system of elements, can be written as: 
[Kk] {w} = {P} OA 


in which {w} now contains all nodal displacements and {P} all nodal 
loads. 
Superposition of the individual element stiffness matrices 
is used to obtain the stiffness matrix for the complete structure. 
Restrained displacements, boundary conditions, may be taken 


into account by greatly magnifying the corresponding diagonal 
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First-Order and Geometric Stiffness Matrices 
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stiffness influence coefficient. A number in the order of 10°° was 
suggested by Aas—Jakobsen. 

The axial force P must be fice in order to evaluate the 
element matrix [K. | in-Hig. 2.9, the axial force as usually not 
known in advance, and an iterative procedure must be used. In the 
first cycle P is chosen equal to zero and the first order forces are 
calculated. In the second cycle the axial forces found in the first 
cycle are used. 

Usually the axial forces are practically not influenced by 
the second-order effects, such that two cycles are generally sufficient. 

In checking the accuracy of the method it was found that the 
calculated moments for frames permitted to sway were in excellent 
agreement with the exact ones when the columns were represented by 
one element. For braced columns and frames the errors in moments were 
as large as 10 percent when the columns were represented by one element 


but vanished when the columns were divided into two elements. 
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CHAPTER 3 


SELECTION AND SECOND-ORDER ANALYSIS OF BUILDINGS 


aoe LoGroduction 

In this chapter eleven building frames are derived. The 
characteristics and loading conditions of each frame are discussed. 
Also, the "exact" method of second-order analysis used to check the 


approximate results is presented. 


3.2 General Information 
3.2.1 Material Characteristics 
All frames studied were assumed to be constructed of normal 
weight concrete with a strength of 4000 psi. The elastic modulus of 


(7) 


the concrete was given by the ACI Code expression, Equation 3.1: 


Eo ay SAO 88) f Sei 


The moment of inertia for all members was based on the 
overall dimensions of the section. To represent the effect of cracking 
at ultimate conditions the moment of inertia of the columns and beams 
was multiplied by 0.8 and 0.4, respectively. These values were arbitrarily 


(9) (14) 


selected based on studies by Hage and the values proposed by Kordina 
The moments of inertia of the shear walls at ultimate were multiplied 
by 0.4 over their entire height, which, in effect, corresponds to the 
assumption that the wall is cracked over its entire height. 

At service loads the members were assumed to be uncracked 
and the gross moment of inertia was used. 


Creep deformations were not considered since the lateral 


loads were all of short duration. 
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3.2.2 Loading Conditions 
All buildings were analyzed under combined dead and live load. 
The dead load values were Pies as: 

i) Frame self-weight. 

Ti) Superimposed dead load, including weight of slabs, mechanical 
equipment and roofing weighing 80 psf on the roof and 100 
psf on the floors. 


Dis) Exterior walls weighing 50 psf of surface area. 


The live load values were taken as: 
1) Rectangular wind distribution of 20 psf. 
Li) Superimposed live load of 30 psf on the roof and 100 psf on 
the floors. 

The frames were analyzed at service and ultimate conditions. 
The loading cases were: 

Service: 1.0(D + L + W) 

UltamatecamOryocliasDo + ln/ ial 7M) 
wiete tne load combination. factor, O02 /5; andthe load factors, 1-4-and 
1.7, for the ultimate load case are those of the 1971 ACI Building Code 
and! GsA7A23 39-1973). 

Wind loading was applied as a concentrated load at each floor 
level applied along the left column line. The right column line is 
therefore denoted as the leeward face. Exterior walls, as noted, also 
formed a dead load component. This loading was lumped into a point 


load applied vertically at each floor level. 


3.3 Details of Frames Studied 
Figures 3.1 to 3.5 present the geometry and member sizes of 


the frames studied. Frames one to five and eight to eleven are assumed 
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to be typical interior frames in buildings containing no bracing not 
already included in the typical frames. In all cases the transverse 
spacing of the frames in the building is assumed to be twenty feet, 
which means the wind pressure, dead and live loads act over a twenty 
foot tributary width. 

Frame type one is a twenty storey, two bay structure as shown 
in Figure 3.1. The sizes of the beam and column cross-sections are 
given immediately to the right of the frame elevation and apply to all 
frames in that figure. The larger of the two column sizes given denotes 
the size of the interior column. The smaller dimensions refer to both 
exterior columns. 

With frame type one it was desired to represent a fairly 
regular structure with a uniform change in stiffness along ies heient.. 

Frame type two, also shown in Figure 3.1, is essentially the 
same as type one. It was desired, however, to investigate the effect 
of a tall and more flexible ground SEOr ee enonee the first storey was 
enlarged from twelve to eighteen feet. 

Figure 3.2 illustrates frame type three which was chosen to 
represent a structure of essentially constant stiffness. One change 
in column sizing was introduced at mid-height of the building because 
it was desired to investigate the effect of the uniform stiffness 
yet represent a feasible structure, 

Frame types four and five are shown in Figure 3.3. Both 
represent shear wall frames with exterior columns and beams the same 
as those of frame type one. The size of the shear wall in type four 
was assumed to be such that its moment of inertia was five times the 


sum of the moments of inertia of the exterior columns. For purposes 
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of analysis the wall was input into the program as a straight line 
member ignoring the effect of the wall width. Also, as previously 
mentioned in section 3.2.1, the mouene of mere of the shear wall 
at ultimate was taken as 0.4 times the gross moment of inertia over 
the total height of the wall. Since some of the wall would probably 
not crack in all storeys, the effect of wall cracking is probably 
overestimated. 

Frame type five was chosen to investigate the effect of a 
heavier wall than that found in type four. In this stiffer structure 
the wall moment of inertia at service loads was taken as fifty times 
the sum of the moments of inertia of the exterior columns. 

Figure 3.4 shows the plan of buildings six and seven. The 
plan shown is one half of a symmetrical structure. The building is a 
flat plate structure with dimensions as shown. The slab thickness is 
taken as seven and one half inches at every floor level. Column and 
shear wall sizing is constant for the total height of the building. 

In modelling this twenty storey building, the frames were 
lumped together and linked to the shear wall-frame as shown in Figure 
3.5. The six inch links were pinned at each end and made axially stiff 
to represent a rigid floor diaphragm. The shear wall was modelled as 
a vertical member with the correct axial stiffness along the center- 
line of the wall and with essentially rigid arms extending from the 
wall centerline to the ends of the floor beams in each floor. This 
allowed the use of a common plane frame analysis program. 

Frame type six is shown in elevation in Figure 3.5. Frame 
type seven was essentially the same as type six only that the shear 


wall was terminated at the seventeenth storey and replaced by two lines 
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of columns above this level. This was done to investigate a frame 
subjected to the whipping action which might develop due to the change 
in stiffness. 

Frame types eight to eleven, shown in Figure 3.1, represent 
a series of special case structures. These four frames are based upon 
the geometry and member sizes of frame type two incorporating various 
modifications to introduce soft storeys. 

Frame type eight was created from type two by deleting the 
second floor, thus creating a thirty foot first storey. Frame type 
nine was formed by deleting the tenth storey of type two. Frame type 
ten was the most radical frame analyzed and was created by deleting 
the second, tenth and top floors of type two. In each of these three 
buildings it was desired to study the effect of the discontinuities. 
It should be noted that these buildings were highly artificial since 
no attempt was made to stiffen the structure in the vicinity of the 
missing floors. 

Frame type eleven was identical to type two in construction. 
However, increased loads, were applied to the roof and top two floors. 
The roof and floors of the nineteenth and twentieth storeys were subjected 
to a superimposed dead load of 480 psf and 500 psf, respectively. 

This additional load served to represent mechanical floors, and tended 


to increase the P-A effect. 


3.4 "Exact" Second-Order Analysis 


To provide a check on the approximate methods of analysis 
it was required to perform an accurate second-order analysis on each 
frame. The negative bracing member method, as discussed in section 2.3.2, 


was used to do this. A negative brace was sized for each storey and 
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inserted with positive slope extending from the bottom of the windward 
exterior column to the top of the leeward exterior column. 
The analyses were carried out using the program PFT (plane 


frames and trusses) described by Beaufait et agit e 
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FIGURE 3.2 ELEVATION OF FRAME TYPE THREE 
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FIGURE 3.3 ELEVATION OF FRAME TYPE FOUR 
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CHAPTER 4 


APPROXIMATE ANALYSIS 


aoe Introduction 

Section 4.2 of this chapter presents the theory behind three 
approximate methods of second-order analysis which will be compared to 
accurate P-A analyses. These are the Fey method, the moment magnifier 
method and the amplified lateral load method. The procedure used to 
apply each method to the frames described in Chapter three is also 
discussed. The errors associated with each method are presented in 
graphical form and is discussed in section 4.3. 

The column moments obtained by each analytical procedure are 


summarized in figures presented in Appendix A. 


4.2 Theory and Application of Analyses 


4.2.1 Fey Method 


From the theory of section 2.2 it is seen that the second- 
order deflections may be estimated from the first-order deflections 


by the equation: 


A 
a 
A, SPA 4.1 
Pte 1 
Hh 
where, 
xP = total vertical load on the storey. 
H = total lateral shear on the storey. 


(18) 


This equation was derived by Poe Ob Goldberg and in a slightly 


C2} 


different form by Parme Equation 4.1 was used to estimate second- 
order moments in the manner described in the rest of this paragraph. 


Initially, a first-order analysis was conducted on each frame at 
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ultimate conditions. The relative lateral floor displacements then 
gave the A, value for each storey. Knowing the values of XP, H, and 
the storey height, h, the value of A, was calculated for each storey 
from Equation 4.1. The A, values for each storey then gave a displaced 
structure which approximated the final shape including the second- 
order effects. Based on these final approximate deflections, column 
end shears were determined as outlined in section 2.3.1. Using these 
end shears, sway forces equal to the difference in end shears taken with 
the appropriate sign were calculated. Finally, these sway forces were 
added to the original lateral load, and the total lateral loads were 
used to estimate the second-order forces, moments and deflections. 
Both the lateral and vertical load moments can be obtained together 
in this analysis if desired. 

Use of the Fey equation (Eqn. 4.1) avoided the necessity 
for several cycles of iteration to compute A, and the second-order moments 
and forces. 

The analysis described above was carried out at the ultimate 
load level. The results from this analysis will be compared to an 


"accurate" second-order analysis in section 4.3. 


4.2.2 Moment Magnifier Method 


Considering the column of Figure 4.1, the bending moment, 


M,, at the base of the column shown by the dashed lines is: 


De 


M, =" Bh PA, ey 


Tf the critical load for the column is denoted by oe and Q = LP/P 


then A, is given as: 
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FIGURE 4.1 


FORCES ON DEFLECTED COLUMN 
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and the second-order bending moment is given as: 


i. tl 
My) = Hh + IPA, Gu5 | 4.4 


From the theory of section 2.2 the value of Q may be approximated by: 
mts 1 4.5 


and the second-order moment, M,, can be approximated by: 


i 


u meee Me 
i eae Ga: 4.6 


Thus, for a given lateral loading pattern leading to first-order frame 


moments, Mi > the total second-order moment, M is given as: 


De? 


where 


1 = ———— 
1-Q DPA, oes 


as ete 


Initially, two first-order analyses were conducted on each 
frame at ultimate conditions, one with vertical frame loading only, the 
second with lateral loading only. From the first-order lateral load 
analysis a value of Q, as given by Equation 4.5, was determined for 
each storey. The moment magnifier, 6, for each storey was then computed 
by Equation 4.8. The column end moments obtained from the first-order 
lateral load analysis were then magnified. These magnified moments 
were added to those from the vertical load analysis to obtain final 


end moments. 
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4.2.3 Amplified Lateral Load Method 

This method attempts to amplify the lateral load such that 
convergence upon the final Lae ca moments is obtained in one 
analysis. The lateral loads to be used in the analysis are taken as 
1/(1-Q ) times the actual lateral loads where o = (2P /H) (A /h) . 
Initially, Qa is based on an assumed drift (A fh) - If the value of 
A fh in any storey computed in the frame analysis exceeds the assumed 
value, the designer would have to rerun the analysis with still larger 
lateral loads, or preferably, he should revise the structural framing 
to reduce A fh to the desired value. 

The serviceability reasons it is necessary to know the total 
(ie..+second=order) deflections at service loads. Reference (9) has 
shown that the EI at service loads is roughly 1.7 times that at ultimate 
loads. For the wind loading case the ACI load factors are 1.05 Dead + 
Lato abive - 1.28. Winds tThus, the.first-order.lateral Load deflection 
at ultimate is 1.28 x 1.7 = 2.2 times the first-order service load 
deflection. However, from Equation 4.1, the relationship between 


service load and witimate load deflections..can_ be estimated as: 


A, Noh 
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where the subscripts s and u refer to service and ultimate loads 


respectively. If we let (A, ,/h) = c(A, /h) » H = 25 Ho 
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A A 
($4) 1(=2) = K and? = G0 tO) 1.28) P then from Equations 4.9 and 
4510 we get: 
ce ds sal Yu Peat 
K Q ° 


ee ee 
R05 etbot 1206) x. 2.7 


For the eleven frames studied the value of K varied Boon 204 
to 2.51, with an average value for frame buildings of 2.11 and for 
buildings with walls of 2.42. The latter value is expected to be on 
the high side in practice since in the analyses the wall was assumed 
cracked in all storeys. In applying this method to the eleven frames, 
the first-order deflection at ultimate was taken as 2.2 times the 
first-order service load deflection (K = 2.2), which was roughly the 
average value for all eleven buildings. Establishing oo =), 25 wile 
is a reasonable upper limit on the allowable value of Q for a storey 
and solving Equation 4.11 for a lower bound, gave c = 0.40. Thus, the 
ratio between the total second-order drift at service and ultimate 
conditions is estimated to be 0.40. 

In carrying out a design, the designer would determine the 
service load sway index required for occupant comfort or to prevent 
non-structural damage, would divide this by 0.4 and get an estimate of 
A fh and use this to compute an and the amplified lateral loads for 


each storey. The value of a is then computed as: 


@ 2. 4.12 


From Equation 4.12 a value of oe can be determined for each storey. 


The amplified lateral load can then be computed as: 
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In design, only one first-order frame analysis would be needed 
and would be carried out using this lateral load. 

In order to apply this method to the coche considered in 
this thesis, a somewhat different procedure was used because the frames 
had not been designed to meet any particular sway index. Rather than 
arbitrarily assuming a value of A,/h at service loads, the maximum 
A, fh value in any storey of each building was used directly. The 
average second-order service load drift was taken as 0.85 of the maximum 
value based on studies of the drifts in frames one, two, three and four. 
The ultimate second-order drift used in comparing this procedure to the 


accurate analyses was then taken as: 


“ou dee x 0.35 = 
h Cc : h 
A 
: Zs 
Se? Lo See 


4.3 Comparisons of the Three Approximate Methods to the 


‘Accurate'' Analysis 


4.3.1 Comparison of Moments 
The relative error in each of the approximate methods with 
respect to the negative bracing solution is herein calculated for each 
frame and presented graphically. The “error function , “e, “was ‘defined 
as the difference in a column end moment between the approximate and 


"exact" value divided by the total storey moment times 100%. That is: 


M - M 
_ approx exact é 
e = x 1007 Hek5 
oti nes + Mow 


The error function is expressed in this way to show the relative importance 
of the error when compared to the total moment governing the design of 


the storey. 


; aA uh 


ah 


on Miaadew esis +1 a 


= 


amend? wide A cilaeios down cue Sib 


nets voter kale howe. est ra \, 


Py aes aa 
eee ee Ea _ -) 


nindeeet ods, abugl ero i ahgh to ou oly ¥ & | “ o a 
; ae het a ‘ a; 
gant i ale hows “ew ighERehbd rhoks dn om va\. 


Puke 


P : 
a bee 


ork ain Gels ~~ cb,{) ea cele vA ann eT brat tere 


: A i As oh } 
. Lo > ar wi rez 
eet bee aend ad sen: onion ar teea _ nodbuse ane, Oa 
. 4 mn 


Ah7 G2 SwIbHser7, i rl? quiketanby ht ‘ber ‘weak sateen sonatas tam 
v as 


a (ee ‘nove ned ney quienes Te 


i: a ek ele 
i? ce a | é 
a. d nt, Se 
“Wah i” je So. 
p afk A ites | 
5 on esi. 4 | oF © 
aate 
1 a) 


d 
| ait? oy abit aah scab iba Gaze! 464 Fo enont xine ti 
a EY, a ea ay Se re Ib, 


Trae Jp BP QO £608 

iy x, : 
iytw sbodtem stimmixoragys add. 36 Apee tk +*paae weselsz oar 
ee 
cone 208 ietiidiciding atoiat at goliulow griongd svl ingot oh? Me 4 iwag 


bert toh) caw - < ok amend tori” att eiiact Aqesy bo'taseexq, bow a 
win os tamtxotygs sao tres dud 2 exepAreaet fae rasa ‘ ob aon te} % oe a 

is ie 
rant I 


jek dott 200, eum omacon yonere hana old fon, awit ay tf ot 


Zia ot is oe Ww 
y ' 


svsiveiay veaiyot aint Bacar “na a 


40. 


A different e can be defined for each end of the column. To condense 
the presentation, only the error in the larger column end moment was 
computed since this moment would govern the design of the column. 

The results are presented in Figures 4.2 to 4.12. The errors are presented 
for the center-line columns in frames one to five and eight to eleven. 
The windward column moments were not considered as they were relatively 
low and would not govern the design. The errors in the leeward side 
column moments were not plotted since, in general, the errors in these 
moments were found te be about one-half those of the center line columns. 
In frames six and seven the errors were calculated in the column line 
immediately to the left of the shear wall. This corresponded approxi- 
mately to a central column line and allowed investigation of the effect 
of the shear wall. In Figures 4.2 to 4.12 the errors corresponding to 
the Fey method are plotted with solid lines, the moment magnifier 

method with dashed lines and the amplified lateral load method with 
broken lines. The horizontal scales change in some of the graphs. 

As shown by the graphs, the approximate methods appear generally 
conservative. In all frames the moments from the Fey method and the 
moment magnifier method were within two percent with few exceptions, 
while the amplified lateral load method was generally considerably more 
conservative. The errors would tend to be higher than expected in 
practice due to the high degree of flexibility in the frames, which 
resulted in higher than normal values of Q = UPA, /Hh and hence higher 
than normal second-order effects. 

The errors in moments were larger and somewhat more erratic 
in frame type five than in the other frames for two reasons. The 


comparisons were carried out for the moments in the shear wall which 
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resisted the majority of the moments hence errors in these moments 
were a larger percentage of the total storey moment. Secondly, there 
were large discontinuities in the wall stiffness and inaccuracies 
appeared at each discontinuity. 

Frame types six and seven showed good results. The error 
values were small and generally conservative. The fact that these 
frames were stiffer than the other frames studied improved the accuracy 
of the amplified load method. 

For the special case structures, types eight to eleven, 
the Fey method was the most accurate. The moment magnifier method gave 
similar results but became inaccurate in the area of the discontinuities, 
in this case the large storeys. The amplified load method was shown 
to be grossly conservative due to the fact that the calculations were 
based on a large storey drift which occurred in the high storeys. 

In a practical structure this error would be expected to be less since 
the columns and beams in such a storey would be stiffened to reduce 


the sway. 


4.3.2 Comparison of Deflections 


Table 4.1 summarizes the roof deflections for each building. 
As seen in the table, each approximate method generally overestimates 
the "exact" tip deflection. At the same time, however, the Fey method 
and the moment magnifier method give a very good approximation for the 
more Jee frames, types one to five, and the frame-shear wall structures, 
types six and seven. With these two analyses the approximation for the 
special case buildings, types eight to eleven, is also good, with 
errors in deflections of generally less than six percent. The amplified 


lateral load method gives a greater degree of conservatism in estimating 
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the roof deflection, consistent with the column moment approximations. 
The error in the amplified lateral load method was approximately ten 
percent for frames one to seven and Peieed as high as fifty percent for 
the special case buildings, types eight to eleven. As pointed out 
earlier this was due, in part, to the manner in which these frames 
were derived with one or more high storeys without increasing the 
stiffness in these regions. 

Table 4.2 illustrates the average second-order drift for 
each building at ultimate conditions. These values are also compared 
to the maximum drift in any storey for that building at ultimate 
conditions. In building types one to seven the maximum storey drift 
ranged from 1.1 to 1.7 times the average and could be taken roughly 
as one and a half times the average. In buildings eight to eleven the 
maximum storey drift was about twice the average. In the derivation 
of the amplified lateral load method in section 4.2.3 this ratio was 
taken as 1/0.85 = 1.18. As a result the amplified lateral loads tend 
to be on the high side. 

The computed maximum storey sways at ultimate ranged from 
1/323 fon frame six to a maximum of 1/76 for frame ten. These values, 
particularly the second value, are higher than would normally be 
expected in reinforced concrete buildings at Tien os Since the 
accuracy of the approximate methods tended to increase as the frames 
became SP ose the inaccuracies observed in the comparisons made will 


tend to be on the high side. 
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TABLE 4.1 Comparison of Calculated Roof Deflections 


Roof Deflections - Inches 


Ultimate 


Frame 
Type vervice 


A, values from negative brace method. 
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TABLE 4.2 Comparison of Maximum and Average Sway Angles 


A, /B Values 
2u 


Frame 
Type Average Max. Storey 
1 .0041 .0056 (6)* 
2 0043 005% (6)* 
3 .0041 .0069 (3)* 
4 . 0036 005 (6)* 
5 .0031 0042 (Ry 
6 .0023 , 0031" 5C/)% 
7 .0023 .0033 (9)* 
8 .0053 SO L245 
9 0051 20122 1) 
10 .0063 50132) (ay * 
11 . 0046 .0060 (6)* 


*Indicates storey in which maximum 


value occurred. 
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FIGURE 4.3 MOMENT ERRORS IN FRAME TYPE TWO 
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FIGURE 4.4 MOMENT ERRORS IN FRAME TYPE THREE 
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CHAPTER 5 


COLUMN DESIGN FOLLOWING SECOND-ORDER ANALYSIS 


J... antroduction 

Figure 5.1 shows columns with and without lateral displace- 
ments of the ends. If translation is prevented, the deflected shape 
is as shown in Figure 5.l(a). The moments Me and M, are the applied 
end moments while Mor and Mop are restraining moments caused by the 
rotations of the end restraints as the column deflects. Horizontal 
forces, H, are present if the end moments are unequal. At mid-height 
there are secondary moments equal to the axial load times the deflection 
shown shaded. To account for the restraining moments Mor and Map in 
the design of this braced column an effective length less than the 
real length is used to compute the lateral deflections. 

If, however, the column is free to sway laterally as shown 


in Figure 5.1(b), the moments M,, and M, must equilibrate not only any 


A 
horizontal load, H, but also a moment PA. The secondary moments in 
this column can be divided into two components, one due to the additional 
horizontal, reaction or sway force, PA/h, necessary to resist the axial 
force in the deformed position, and the second equal to the axial load 
times the deflection from the chord line, shown shaded in Figure 5.1(b). 
Traditionally these have both been accounted for in design by using 
the elastic effective length factors for the unbraced case in designing 
the column. 

On the other hand, if a second-order structural analysis 


is carried out including the effects of both the applied loads and the 


sway forces, the latter have been accounted for in the analysis and 
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need not be considered a second time in evaluating the effective length. 
Since the maximum moment theoretically may occur away from the ends of 
the column it may be necessary to use a moment magnifier calculation 

to estimate this moment. A procedure for determining whether such a 
magnification is required is presented in section 5.2, and the necessary 
values of the effective length factor, K, and the equivalent moment 
factor; C. are discussed in section 5.3. The proposed column design 
procedure is compared in section 5.4 to the second-order column moments 
in a one storey frame and to 10 columns selected from the 11 frames 


analyzed in Chapter 4. 


5.2 Test of Whether Maximum Column Moment Occurs at End of Column 


(6) 


Galambos shows that the maximum moment in an elastic 


beam-column loaded with an axial load and end moments M, and M is: 


MAX = OM Sei 


ae (Mm, /M)? = 2(M, /M.) cosa 


sind 
where M, = larger end moment in the column, always positive. 
M = smaller end moment in the column, positive if column 


is bent in single curvature. 


a = hv P/EIL = Tv P/P., 


2 
fe Lids 
Po = Euler load = wines 
(¥1)) . ; 
MacGregor has used these equations to derive a test for 


determining whether a given column can be designed for the maximum end 


moment , M.> or whether the moments between the ends of the column will 
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exceed those at the ends. If it is assumed that stability effects can 
be disregarded if Max is) not more: than 1705°M.. Equation 5.2 can be 
a 
solved to determine the combinations ‘of MM, and @ corresponding to 
6 = 1.05. These are plotted with the solid line for 6 = 1.05 in 
Figure 5.2. Columns represented by combinations of M/M, and a? 
falling below this line can be designed for the second-order end moments 
without a further moment magnification. This line can be approximated 
by the equation shown in Figure 5.2. Thus if: 
me Poe 
u 


es ee = 
M ol 3EI1 
a 


) 52 


The maximum moment will always be less than 1.05 times that at the end 


of the column. 


5.3 Derivation of K-Factor 

If Equation 5.3 shows that the maximum column moment occurs 
away from the end of the column, column Jeon should be based on 
amplified moments based on the ACI moment magnifier with the equivalent 
moment factor, Ce taken for the braced frame case using the ratio of 
end moments obtained from the second-order analysis, and with the 
effective length factor K = 1.0. This can be demonstrated by setting 
the moment magnifier from Equation 5.2 equal to the ACI moment magnifier 


given as: 


4 


where, Ph is the ultimate load on the column and F ot the critical load 


equal to m7EI/(KL)*. 
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LIMITS ON COLUMNS WITH MAXIMUM MOMENT AT END OF COLUMN 


60. 


Oke 


and 


Ce = 0.674-0.4 vl De 


Equating Equations 5.2 and 5.4 gives: 


C ; {1+ Ce = 2(M, /M_) cosa 


sind 


which reduces to: 


ES ee eee ee ae 
ag ee = 57 
> ele (M,/M,)? - 204 /M,) cose 


Equation 5.7 was solved for values of M/M,, from -1 to +1 and P/P., 

from 0.03 to 0.63. The values of K computed using Equation 5.7 ranged 
from 0.999 to 1.265. The values of K greater than 1.10 corresponded 

to cases where Ch is underestimated by the ACI equation (columns with 

Q and M/M both low. See Reference 205) If these values were excluded, 
the average K value was equal to 1.05. The absence of any K values 

less than 1.0 indicates that the braced frame effective lengths should 
not be used if the column end moments are known from a second-order 


analysis. 


5.4 Verification of Design Procedure 


To verify the column design procedure, presented in sections 
5.2 and 5.35. the mrame of ee 5.3 was analyzed with a negative 
bracing member. The loading was as shown and the column under 
consideration was broken into four equal length segments to allow 


calculation of the deflections and moments at points between the ends 


of the column. 
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The frame was assumed to be constructed of normal weight 
concrete with a strength of 4000 psi and elastic modulus as given by 
Equation 3.1. The column and beam moments of inertia were 15617 in’ 
and 16384 in’, respectively. The negative brace area was given by 
Equation 2.17. 

Initially, the frame was analyzed as shown in Figure 5.3 
to obtain deflections and end moments in the windward column. From 
the resulting deflections it was possible to calculate the relative 
deflection from a chord line joining the upper and lower ends of the 
column. The frame was then reanalyzed under the same loading but with 
the node points in the position resulting from the first cycle. Three 
cycles were required to obtain convergence. The moments in the windward 
column of the frame are shown in Figure 5.4. . During the iterations the 
moments at mid-height of the column increased and those at the ends 
decreased. The decrease in the end moments occurred because the lateral 
deflections caused a slight decrease in the stiffness of the column. 

Equation 5.3 may now be used to test whether design of the 
column may be based on the larger end moment or whether a magnification 


is required. The ratio of column end moments is found to be: 


The right-hand-side of Equation 5.3 is: 
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Hence, the test shows that the maximum moment occurs between the ends 
of the column and design must be based on an amplified moment. 

Design of this column by the ACI code would require calculating 
Max as given by Equation 5.1. 

For the frame of Figure 5.3, based on the second-order end 


moments from the third cycle of iteration, Ch would be: 


i BBR 3 
Ch = (Ose 0.4057 


) f= 5.160 


The critical buckling load assuming an effective length factor equal 


Lo) ..0 would. be: 


(3.14)7 x 3605 ksi x 15617 in* 


cr CUSOrx 144° in) AOU SENS. 


Solving Equation 5.4 for the moment magnifier gives: 


eee ly OS es 
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26769K 


Finally, the column design moment would be given by Equation 5.1 as: 
MAX =e 1.23) x 15. 45 ft-Kki = 216.9 stk 


This compares to a maximum column moment of 19.7 ft-K from the iterative 
analysis. Thus, the ACI procedure with K = 1.0 underestimated the 
maximum column moment by 3.8 per cent. Had the braced frame effective 
length of 0.81 been used in the calculations, the magnified moment 
would be estimated as 17.6 ft-K or 10.9 per cent low. This suggests 
that the hinged end effective length (K = 1.0) gives a better estimate 
of the magnification than the braced effective length. 

A number of columns were chosen from the eleven frames 


studied in Chapter four to further check this design method. Table 5.1 
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gives the second-order end moments of ten single curvature columns and 
the ultimate vertical load on the storey where that column is located. 
In calculating the second-order ene: along the columns the total 
vertical storey load was artificially assumed to act on the column in 
question to represent an extreme hypothetical condition for which the 
second-order frame analysis was applicable. As shown in Table 5.1, 
Equation 5.3 indicated that six of the ten columns required additional 
moment magnification to estimate the maximum moments in the columns. 

The moments computed using a moment magnifier analysis are 
compared in Table 5.2 to those from the second-order analysis. In 
both calculations the entire axial load was assumed to act on the 
column in question. It can be seen that columns one, five, six and 
eight, which Equation 5.3 identified as not requiring additional moment 
magnification, indeed did not need moment magnification. All the 
columns identified as requiring moment magnification, did need magni- 
fication except column ten which was marginal but did not require 
magnification. 

The values of the maximum column moment in Table 5.2 computed 
using K = 1.0 were conservative in all cases while those based on 
the braced frame K were unconservative in all cases when compared 
to the column moments computed in one cycle of iteration by the 
procedure described earlier. This supports the observation that 
if the column end moments are known from a second-order analysis, 
the effective length factor should be taken as K = 1.0 and the 


equivalent moment factor, Cw should be based on the braced frame case. 
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TABLE 5.2 Comparison of Moments Computed Using Moment Magnifier 


to Second-Order Moments 


"MAX 


Actual 


* Maximum moment at end of column. 
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CHAPTER 6 
SUMMARY AND CONCLUSIONS 


In this study a column design method, based on a second- 
order analysis, was proposed. It was shown that the maximum moment 
in a single curvature column could be estimated by applying a moment 
magnifier to the larger end moment. 

Results indicated that the effective length used in 
calculating the moment magnifier be taken equal to the actual 
column length. The use of the braced frame effective length gave 
approximations to the maximum column moment which were unconservative. 

Three approximate methods of second-order analysis were 
presented and applied to eleven different building frames. 

The Fey method and the moment magnifier method were shown 
to give essentially acceptable approximations to an "exact" second- 
order analysis. The amplified lateral load method, however, was 
grossly conservative in some of the cases studied, but became comparable 
to the Fey method and moment magnifier method for the stiffer frames. 
Difficulties also arose in this method if the sway deflection angles 


varied significantly from storey to storey. 
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APPENDIX A 


COLUMN END MOMENTS FOR FRAMES STUDIED 


2a we LntLoduction 

In this Appendix are tabulated all column end moments for 
each frame studied in this thesis. In each figure the upper right 
hand box attached to the column indicates the top end moment. The 
bottom end moment is listed in the lower left hand box. The first 
value in each box represents the "exact" second-order end moment. 
The second, third and fourth entries give the results of the Fey 
method, the moment magnifier method and the amplified lateral load 
method, respectively. The arrow associated with each box indicates 


the direction of the end moment on the column. 
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FIGURE A.3 FRAME TYPE THREE 
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FIGURE A.4 FRAME TYPE FOUR 
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FIGURE A.5 FRAME TYPE FIVE 
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FIGURE A.6(a) FRAME TYPE SIX 
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FIGURE A.6(b) 
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FIGURE A.7 (a) FRAME TYPE SEVEN 
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FIGURE A.10 FRAME TYPE TEN 
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